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OF GENERAL TORAL ENDOMORPHISMS 

MICHAEL BAAKE, EIKE LAU, AND VYTAUTAS PASKUNAS 

Abstract. It is well-known that the Artin-Mazur dynamical zeta function of a hyperbolic 
, or quasi-hyperbolic toral automorphism is a rational function, which can be calculated in 

terms of the eigenvalues of the corresponding integer matrix. We give an elementary proof of 
this fact that extends to the case of general toral endomorphisms without change. The result 
is a closed formula that can be calculated by integer arithmetic only. We also address the 
£^ functional equation and the relation between the Artin-Mazur and Lefschetz zeta functions. 
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1. Introduction 

Q 

Any (i-dimensional toral endomorphism is represented by an integer matrix, M G M&t(d, Z), 
with action mod 1 on the d-torus T d ~ M^/Z^; see [2] and \15\ Ch. 1.8] for background and 
[H Ex. 1.16] for an illustration. Important aspects of the dynamical system (T d ,M) are 
related to its periodic orbits and their distribution over T d ; compare |11[ I24j. The Artin- 
Mazur [3j dynamical zeta function provides a generating function for the orbit counts that is 
interesting both from an arithmetic and from a topological point of view |13t [20] . The latter 
was also Smale's approach [22], who related the Artin-Mazur and Lefschetz zeta functions of 
a hyperbolic toral automorphism and calculated both in terms of eigenvalues. 

In this note, we explain a different approach via elementary geometry and linear algebra, 
which bypasses more advanced topological methods as well as the need to calculate eigenval- 
ues. A key observation is that all arguments apply to general toral endomorphisms without 
additional effort. We also treat the connection between the Artin-Mazur and the Lefschetz 
zeta function and their functional equations. As we have learned along the way, most argu- 
ments we use appear already in the literature, notably in [13], but at least their combination 
seems to be new. Also, we make several steps explicit to facilitate their computational use. 

For M £ Mat(<i, Z) and m > 1, let a m be the number of isolated fixed points in T d of the 
m-th iterate M m . The starting point of our considerations is the identity 

(1) a m = |det(l-M m )|. 

This formula is well-known |24[ 3] when no eigenvalue of M m is 1; it then follows from 
counting the number of points of Z rf in a fundamental domain of the lattice (1 — M m )Z rf . 
Otherwise, (pQ) is true because both sides are zero. Indeed, since the fixed points of M m form 
a closed subgroup of T d , they are either all isolated, or they form entire subtori of positive 
dimension; see [H Appendix] for a detailed discussion of the subtorus case. Incidentally, when 
no eigenvalue of M m is 1, a m is also the Reidemeister number of a toral endomorphism, see 
[T3| Thm. 22 and p. 33], while, in general, a m is its Nielsen number. 
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Following [3], the Artin-Mazur zeta function of a general M € Mat(d, Z) is defined as 

oo oo 

(2) Q M (z) : = exp(£ ^z m ) = ]J (1 - z m )~^ . 



m 

-1 m=l 



Here, the exponents c m of the Euler product representation are well-defined integers, see 
Proposition [2] below. An explicit representation of Cm( z ) as a ra ti° na l function is given below 
in Theorem [TJ We stress that at least for hyperbolic or quasihyperbolic toral endomorphisms, 
this result is well-known by [22] or [13]; our focus is the elementary method. 

A matrix M is called hyperbolic when it has no eigenvalue on the unit circle § . Such 
toral automorphisms are expansive [24, p. 143]. Note that M may possess eigenvalues on 
S 1 other than roots of unity (for instance, if one eigenvalue of M is a Salem number; see 
|23|. [5] for examples). Integer matrices without roots of unity in their spectrum constitute 
the quasihyperbolic cases, compare |23] and references therein, where formula ([1]) still counts 
all fixed points. For quasihyperbolic matrices M, the exponents c m are the cycle numbers, 
which are related to the fixed point counts via 

(3) a m = s £ j lc l and c m = — J^/x(y) a t . 

£\m t\m 

This follows from a standard application of Mobius inversion; compare [191 16]. 

When roots of unity are among the eigenvalues of M, the Euler product still exists, with 
the same relation between the counts of (isolated) fixed points and the exponents c m , though 
the latter can now be negative. Let us briefly illustrate this phenomenon in one dimension. 
Endomorphisms of T 1 ~ S 1 are represented by multiplication (mod 1) with an integer n. The 
dynamical zeta function reads Co( z ) = V(l — z ) an< ^ 

(4) CM = 

1 — \n\z 

for n^O, due to our Theorem [T] below (or a simple direct calculation). For n = — 1, we get 
C-i(z) = (1 — z 2 )/(l — z) 2 , thus c 1 = 2 and c 2 = —1, while c m = for all m > 3. The negative 
c 2 corresponds to the fact that the two isolated fixed points of the map fail to be isolated for 
any even iterate, which is the identity. 

Finally, let us note that our arguments extend to the case of nilmanifolds X = G/T 
considered in |13} Sec. 2.6], where G is a simply connected nilpotent Lie group and r a 
discrete subgroup such that X is compact. Namely, any endomorphism ip of r extends to an 
endomorphism (p of X, the isolated fixed points of which are counted by | det(l — Lie(<^))|, 
analogously to ([1]). 

2. A RELATED ZETA FUNCTION 

Let us start with the numbers a m := det(l — M m ), which can be viewed as signed fixed 
point counts, and the corresponding zeta function 

oo ~ 

(5) S*(*)=e*p(£^* m ). 

m=l 

In Section [U we will see that this is actually a Lefschetz zeta function, see Eq. (fl6j) below. 
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For A E Mat(d,R), let A k (A) be the induced linear map on the exterior power A fc (M d ). In 
terms of the standard basis of that space, A k (A) is represented by the matrix of all minors of 
A of order k; see Ch. 1.4] for details. This is an integer matrix of dimension f.), with 
A°{A) = 1, A 1 (^4) = A and A d {A) = det(A). 



fe+i 



Proposition 1. For M G Mat(d,Z), we have ( M (z) = TlLo^ 1 - zA k (M)) { 1] . 

Since all A k {M) are integer matrices, Cm( z ) * s a rational function with numerator and 
denominator in Z[z] . It can be calculated by integer arithmetic alone (many algebraic program 
packages have the matrices of minors of arbitrary order k as built-in functions). Also, since 
is never a root of the denominator, the series © for Q M converges uniformly on sufficiently 
small disks around 0. 

Proof of Proposition^ This is analogous to |X3|, Lemma 17]: The assertion is immediate from 
the well-known formula in linear algebra 

d 

(6) det(l-A) = £(-l)*tr(A*(il)) 

together with the power series identity 

tr(A m ) n 1 



(7) exp (J2 



m 

m=l 



det(l - zA) ' 



which is omnipresent in connection with zeta functions of any kind (in particular, it appears 
in the calculation of dynamical zeta functions of shifts of finite type, see [3 |2Q]). We recall 
that (|6|) is proved by evaluating the characteristic polynomial of M at 1, while ((7|) is a simple 
consequence of the relation det(exp(C)) = exp(tr(C)) for square matrices C, together with 
the Taylor series for — log(l — z), which is the case d = 1 of ([7]). □ 

Let us also note that, in terms of the d eigenvalues X 1 , . . . , X d of M, one has the relation 
det(l - zA k {M)) = P k (z) with the polynomials P (z) = 1 - z and 

p k i*)= n (i-*v----\) 

i<i 1 <e 2 <---<e k <d 

for 1 < k < d. This version is useful for the derivation of the functional equation of ( M . 

Lemma 1. If M £ Mat(d,Z) with D := det(M) + 0, one has ( M (l/Dz) = B (( M {z)) { ' 1)d , 
where B = D for d = 1 and B = 1 otherwise. 

Proof. First, a direct calculation shows that 

i x i 



where /3 = 1 and (3 k = D^- 1 ' for 1 < k < d. Note that each prefactor j3 k involves products 
of eigenvalues, but is symmetric in them and thus simplifies to a power of the determinant. 

Next, recall the binomial formula YH=o 0>)( — 1) = ^™.o f° r n > 0, and insert the previous 
polynomial identities into the product expression of Proposition [TJ Our claim follows, because 
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the prefactor that contains z disappears by an application of the binomial formula, while 
the prefactor with the determinants simplifies to the factor B by an analogous calculation; 
compare |13t Lemma 19] and its proof for a related argument. □ 

The special situation for d = 1 is also immediate from C, n (z) = (1 — nz)/(l — z), as the 
determinant is n; compare the example in the introduction. 

3. The Artin-Mazur zeta function 

To derive a formula for the dynamical zeta function, we observe that a m = a m sgn(a m ). 
Hence the signs of all nonzero a m need to be determined. When M is quasihyperbolic, this 
is done in [23] Lemma 2.1], see also the proof of [13, Lemma 15], but the argument works for 
general M G Mat(d, Z), too: We employ the formula 

d 

(8) S m = det(l-M m ) = {](!- A™) 

i=i 

with the A ■ as above. It is clear that neither complex eigenvalues play a role (as they come in 
complex conjugate pairs, and (1 — A m )(l — A m ) = |1 — A m | 2 > 0), nor do eigenvalues A G [—1, 1] 
(because then 1 — A m > 0). The remaining eigenvalues (evs) matter, and one finds 

(9) a m = a m ((-1) # rcal evs <- 1 ) m (-l) # real evs outsidc I" 1 ' 1 ! =: a m 5 m e. 
Inserting this into Cm( z ) an< ^ comparing with Cm( z ) & ves 

(10) Cm(z) = (Cm(Sz)) £ 

for the dynamical zeta function of M. Though ([9]) involves the eigenvalues of M, the signs 5 
and e can once again be obtained by integer arithmetic alone. When no eigenvalue of M is 
±1, they are simply given by 

(11) 5 = sgn(det(l + M)) and e = <5sgn(det(l — M)). 
In general, the signs can be defined by the one-sided limits 

5 = lim sgn(det((l + a)t + M)) and e = 5 lim sgn(det((l + a) I - M)), 

which can be evaluated explicitly as follows. Factorise det(xt — M) = (x — l)°~(x + l) T Q(x) 
with Q € Z[x] and Q(±l) ^ 0, where the non-negative integers a, r are unique. This implies 
det(xl + M) = (x — l) T R(x) with R E 7L\x\ and -R(l) ^ 0. Consequently, one has 

'det(xl - M) 



( 'det(xl + M) 
(12) o = sgn 1 



x=l 



and s = 5 sgn 



(x - iy 

which is used in our sample program in the appendix. 
Let us summarise the result of our derivation so far. 



(x - l) a 



x=l 



Theorem 1. Consider a general toral endomorphism, represented by a matrix M € Mat(cZ, Z). 
The associated Artin-Mazur zeta function, defined in terms of isolated fixed points, satisfies 

d 



CAf(^)=n det ( i - fcAA; ( M )) £ 



k=0 
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where the signs 5 and e are given by Eq. (jlip when ±1 is not an eigenvalue of M , and by 
Eq. (|12|) in general. In particular, Cm( z ) ^ s a rational function. When no eigenvalue of M is 
a root of unity, all fixed points are covered this way. □ 

For quasihyperbolic toral endomorphisms, this result follows from an analogous result for 
nilmanifolds [13\ Thm. 45], which is proved by Reidemeister-Nielsen fixed point theory, while 
the case of hyperbolic toral automorphisms is already treated in [22\ Prop. 4.5] in a slightly 
different formulation. Theorem [1] covers the special cases of automorphisms for d = 2 from 
|15l [TT] IS]. Let us note that, in the hyperbolic case, the rationality of Q M can also be seen as 
a consequence of the general rationality result [T7] proved by Markov partitions. Of related 
interest is the approach of [18], which connects the problem to an interesting class of % d - 
actions. 

Since 5 and e are signs, in particular 5 = 1/5, the functional equation for £ M is now 
immediate from Lemma Q] and Theorem [TJ 

Corollary 1. When D := det(M) / 0, one has ( M (l/Dz) = B e (( M (z)) ( ~ 1)d with e from 
Eq. (|12p . where B = D for d = 1 and B = 1 otherwise. □ 



Let us also mention that 

zC'm(z) 
Cm(z) 



(13) a - 



z m 



is the ordinary power series generating function of the sequence (a m ) meN , which is still a 
rational function. Its radius of convergence g M is always positive (it is the absolute value 
of the smallest root of the denominator of (|13p in reduced form). Thus, when linim^oo m+1 
exists, 1/ ' g M is the asymptotic growth rate of the fixed point counts, which provides a simple 
alternative to the approach in [23] . The limit exists precisely for hyperbolic (and hence 
expansive) endomorphisms, as follows from [101 Thm. 6.3]. The ratio as a growth measure is 
also employed in [16, Thm. 16], where the case of unimodular roots is briefly discussed, too. 

Finally, we observe that the Artin-Mazur zeta function of a general toral endomorphism 
can be written as an Euler product. 

Proposition 2. For any matrix M £ Mat(<i, Z), the associated Artin-Mazur zeta function 
Cm( z ) has an Euler product representation ([2]) with uniquely determined integers c n . 

Proof. The Euler product representation ([2]) is equivalent to the relations ([3]). These define 
rational numbers c n which we must show to be integers. When M is quasihyperbolic, this is 
true by their geometric interpretation as cycle numbers. The general case (including d = 1, 
which also follows from (J3D) can be proved by a deformation argument as follows. 

Fix n and recall that c n is linear in the a^, while = e> fc edet(l — M k ) from with the 
signs 5 and e. Hence, c n is a polynomial with rational coefficients in the entries of M, which 
can be written as c n = eP(M). Here, P itself depends on 5, but neither on e nor on M. For 
another matrix M' , let 5' , e' and c' n denote the associated signs and numbers. When 5' = 5, 
we thus have c' n = e'P(M'). 

Let v > 1 be a common denominator of all coefficients of P and define the diagonal dxd- 
matrix N = diag(2£, 3, 4, . . .), which shares the sign 5 with M. Consider now M' = M+u r N = 
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v r (N + v~ r M). Since N has distinct real eigenvalues, the eigenvalues of N + u~ r M, for 
sufficiently large r, are real and close to those of N. Then, M' is hyperbolic with 8' = 6 by 
construction; in particular, c' n = e'P(M') is integral. Since the difference P(M) — P(M') is 
integral as soon as r > 1, it follows that c n = eP{M) is integral, too. □ 

4. Interpretation as Lefschetz zeta function 

Suppose that X is a compact differentiable manifold, assumed orientable for simplicity, and 
/ : X — ► X is some differentiable map. In this situation (and also more generally), there is 
a fixed point index /jGZ which satisfies the Lefschetz trace formula 

(14) i f = J2(- 1 ) k Hf*\H k (x;®)), 

k 

see \12\ Prop. VII. 6. 6] or [8j Thm. 12.9]. Here, Hk(X;Q) denotes singular homology with 
coefficients in Q. It is a finite-dimensional Q-vector space in our situation, on which / acts 
by functoriality. 

When all fixed points of / are isolated, we have If = J2xeFix(f) V( x )' where if(x) G Z is 
the local index of / at x. If x is a regular fixed point, meaning that 1 is not an eigenvalue of 
the tangential map T x (f), the local index is given by 

i f (x) = sgn(det(l - T x (f))) e {±1}. 

The Lefschetz zeta function associated to / can be defined as 

C/^)=ex P (E-^ 
n>l 

This definition seems to appear first in [22] ; see also [13] . By using the identity (|7|) , the trace 
formula (I14h applied to all iterates of / implies that Cf( z ) is a rational function, 

(15) C/(^)=n det ( 1 -^*l^(^;Q)) ( " 1)fc+1 - 

k 

Let us now assume that X = T d = M. d /7, d as above, and that / is given by an arbitrary 
M G Mat(d, Z). In this case, with the zeta function C, M of Proposition [lj we have 

(16) ( M (z)=(t(z), 

including a correspondence of all related formulas (for the case of hyperbolic toral automor- 
phisms, this was noted in [131 P- 86, lines 7 and 16]). Let us sketch a possible line of argument. 
First, it is well-known that the Kiinneth formula {8, Thm. 3.2] gives an isomorphism 

(17) H k (X;Q) A k (Q d ) 

such that the action of / on corresponds to A fc (M). This identifies the right hand sides of 
© and (fI3j) , and similarly for Proposition [1] and Eq. (f!5|) . It follows that the left hand sides 
of the corresponding pairs of equations are equal as well, that is 

(18) // = det(l - M), 

and similarly Ifn = det(l — M n ); this appears also in [22j Prop. 4.15] and in [9]. 
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A direct proof of (|18p without using the Lefschetz trace formula can be done as follows. 
Assume first that 1 is not an eigenvalue of M. All fixed points x of / are then regular, with 
the same local index ij(x) = sgn(det(l — M)), because the tangent space T X (X) can be 
identified with M. d , where the action of / is given by M. Thus (|18p is immediate. When no 
eigenvalue of M is a root of unity, the same applies to all iterates of /. 

For arbitrary M, we may use the following lemma. 

Lemma 2. Assume X to be a compact Lie group of dimension d and f : X — > X a dif- 
ferentiable map. Let g{x) = x ■ f(x) . Then, g acts on the 1-dimensional Q-vector space 
Hd(X;Q) by the scalar If. 

Granting the lemma, (|18p follows easily because, in the torus case, g is given by 1 — M, 
and we have Hd(X; Q) = A d (Q d ). We only sketch the proof of Lemma[2]and leave the details 
to the reader. The fixed point index can be defined as the homology intersection product in 
XxX of the graph of / and the diagonal, If = [Fx] ■ [A]. Since the automorphism of X x X 
given by (x, y) i— ► (x, x-y~ l ) acts on the orientation by (— l) d , we get If = (—l) d [T g ] • [X x {1}]. 
The assertion follows by a straightforward computation based on decomposing [T g ] according 
to the Kiinneth formula for XxX. 



Appendix: A sample program for calculating Cm 

One can implement the explicit zeta function formulas of Proposition [1] and Theorem Q] in 
a simple Mathematica® program as follows. 

Clear[tilzeta, zeta, ord,sig,tau, del, eps, dim, one]; 
tilzeta[mat_]:=(dim=Length[mat]; Factor[Product[Det[ 

ldentityMatrix[Binomial[dim,k]] - z Minors[mat,k]]'~((-l)' s (k+l)),{k,0,dim}]]); 
ord[pol_, x_]:=(tmp=pol; i=0; While[(tmp /. z^x)==0, (i ++; tmp=D[tmp,z])]; i); 
zeta[mat_]:=(one=ldentityMatrix[Length[mat]]; pol=Det[z one — mat]; 

sig=ord[pol,l]; tau=ord[pol, — 1]; 

del=Sign[Factor[Det[z one + mat]/(z— l)"tau] /. z— >1]; 
eps=del Sign[Factor[pol/(z— l)"sig] /. z— >1]; 
Factor[(tilzeta[mat] /. z— > del z)"eps]); 

The input is an integer matrix, in the standard format of a double list. The calculation is 
exact and reasonably fast for small dimensions, and can be used up to dimension 8 or 10 say. 
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